Origins of low- and high-pressure discontinuities of T c in niobium 
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The discontinuities of T c in Niobium under pressure are examined by means of the pseudopoten- 
tial plane-wave implementation of the electron-phonon coupling calculated from density-functional 
perturbation theory. Both low- and high-pressure discontinuities of T c have their origin in the Kohn 
anomalies and are caused by the low-frequency phonons, but the mechanism leading to the disconti- 
nuities is different in the two cases. The low-pressure anomaly is associated with a global decrease of 
the nesting factor in the whole Brillouin Zone and not to a visible change in the band structure. The 
high-pressure anomaly is instead connected with a well-pronounced change in the band structure. 

PACS numbers: 71.15.Mb, 71.18.+y, 63.20.Kr, 74.62. Fj 



I. INTRODUCTION 

Niobium is a superconductor with a quite high crit- 
ical temperature, T c — 9.25 K, for a simple metal. The 
experiments under pressure by Struzhkin et a/i show dis- 
continuities of T c at about 5 GPa and at 50-60 GPa. The 
low-pressure discontinuity manifests itself as an increase 
of T c by about 1 K. The high- pressure anomaly is associ- 
ated with a decrease of the critical temperature. To date, 
the nature of these pressure- induced discontinuities is not 
clear. Previous theoretical studies^^^ agree in attribut- 
ing the high-pressure anomalies to some visible change in 
the band structure. The low-pressure discontinuity of T c , 
however, remains mysterious. The goal of this work is to 
give more information about the nature and the origin of 
the anomalous behavior of T c in niobium under pressure. 
In particular, we want to understand the role of Kohn 
anomalies 6 of the phonon spectra. Kohn anomalies are 
known to drastically change the critical temperature in 
superconductors and are believed to play a role in all 
body-centered cubic (bcc) metals Nb, Mo, V, Ta. 7 

Therefore, we look to the details of the electronic struc- 
ture and dynamical properties of Nb at eight pressures 
in the range from -16 GPa to 78 GPa. We show that 
Kohn anomalies are responsible for both discontinuities, 
but the origin of the low-pressure Kohn anomaly is dif- 
ferent from that of the high-pressure one. Both can be 
identified by a closer study of the Fermi-surface nesting 
and of the band structure. 

The accurate calculation of the electron-phonon cou- 
pling A and of the spectral function a 2 F is crucial for our 
problem. To this end we use density-functional perturba- 
tion theorj&S (DFPT) in a pseudopotential plane-wave 
approach. The Eliashberg function and the nesting fac- 
tor require an integration of the double delta over the 
Fermi surface, which needs to be done with a high nu- 



merical accuracy. We give a few advices for an efficient 
calculation of the electron-phonon coupling. Some of the 
technical details, however, can be used in general for cal- 
culations of other properties which require an accurate 
numerical integration with the delta function. 

This paper is organized as follows: In the next Section 
we remind the physical definitions and give some details 
of the calculation of electron-phonon interaction coeffi- 
cients using Vanderbilt's ultrasoft pseudopotentialsri - In 
Sec. Ill, we give the technical details (Subsec. A) and 
present results for several properties under pressure: the 
lattice constant and bulk modulus (Subsec. B), the band 
structure and Fermi surface (Subsec. C), the phonon fre- 
quencies and linewidths (Subsec. D), and the Eliashberg 
function and electron-phonon coupling constant (Subsec. 
E). In Sec. IV, we discuss the origin of the anomalies, 
and we summarize in Sec. V. In the Appendix, we give 
numerical details for the calculation of the Eliashberg 
function. 



II. ELECTRON-PHONON COUPLING 

A. Definitions 

The Hamiltonian for the electron-phonon interaction 
is given in second quantization by 



H, 



el—ph 



kqi/ 



qv,mn Jm n /, t , U \ 
yk+q,k C k+q C k \ u -qv ^ U W ) 



(i) 



where cj^™ and c£ are the creation and the annihilation 
operators for the quasiparticles with energies £k+q.m and 
£~k,n in bands m and n with wavevectors k + q and k, re- 
spectively; and bq U are the creation and the annihila- 
tion operators for phonons with energy w ql/ and wavevec- 
tor q; the matrix element 3k+q™k describes the electron- 
phonon coupling. The coupling constants ffjj+q k define 
the spectral function, a 2 F(cu), its first reciprocal momen- 
tum, A, and the superconducting electron-phonon cou- 
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pling constant, Aq^, by the following set of equations: 
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The quantity N(e F ) is the density of states at the Fermi 
energy, e F , per both spins. 

We introduce, after Allenjii the phonon lincwidth 7 q; ,: 
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which enters the Eliashberg function, a 2 F, and the 
electron-phonon coupling constant, A qi/ , as follows: 
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B. Matrix elements of the electron-phonon 
interactions 

Within DFPT^ the electron-phonon matrix elements 
can be obtained from the first-order derivative of the 
self-consistent Kohn-ShamiS (KS) potential, Vks, with 
respect to atomic displacements, u s r for the s— th atom 
in lattice position R, as: 



qi/,mn ( ^ 



1/2 



<Vk+ q , m |A^|</>k,„>, (8) 



where tpk,n is the n— th valence KS orbital of wavevector 
k and 
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is the self-consistent first order variation of the KS po- 
tential, N is the number of cells in the crystal, and uf J 
is the displacement pattern for phonon mode v^ v : 
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The latter is obtained from the diagonalization of the 
dynamical matrix, <?" s ,(q): 



(ii) 



M s is the mass of atom s, a,/3 denote cartesian coordi- 
nates. 



C. Matrix elements with ultrasoft pseudopotentials 

The use of ultrasoft (US) pseudopotentials (PPs) 10 al- 
lows in many cases a significant reduction of the needed 
plane-wave kinetic energy cutoff, as compared with stan- 
dard norm-conserving pseudopotentials^^^ This en- 
ables a more efficient calculation, at the price of introduc- 
ing additional terms originating from the augmentation 
charges employed in this schemed A detailed description 
of DFPT with US PPs has been given elsewhere by Dal 
Corso^i Here, we only briefly describe terms appearing 
in the electron-phonon coupling. 

With US PPs the KS orbitals, Vk,n> satisfy a general- 
ized eigenvalue problem 
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where the overlap matrix, S, is given by 



S , (ri,r 2 ) = <5(ri-r 2 )+ } , q nm 
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The charge correction q nm , in the above formula, is de- 
fined with the augmentation functions, Q„ m (r — R s ), as 
follows 
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and the projector functions /3„(r — R s ) are specific for 
the type of atom at the position R s , and are obtained 
from the atomic calculations. 

The valence charge density is then computed as 



k 



smn 



(r-R s )(^ k ,i | V>k,i>) 

J] / / d 3 r 1 d 3 r2^ k , i (r 1 )X(r;r 1 ,r 2 )^M(r 2 ), (15) 



where the sum over k and i runs on occupied KS orbitals, 
the kernel 

X(r;ri,r 2 ) = 5(r - n)5(r - r 2 ) + ^ Q nm (r - R s ) 

snm 

xA,(n-R.)/^(r a -R.) (16) 

has been introduced for later convenience. 

The KS selfconsistent potential in the US-PP scheme 
reads 

Vft-s(ri,r 2 ) = V NL {r 1 ,r 2 ) 

+ J d 3 rU e// (r)^(r; ri ,r 2 ). (17) 
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The effective potential, V e ff, contains the local, the 
Hartree and the exchange-correlation (xc) terms 



V eff (r) = V loc (r) + / d 3 n 



p( r i) 
|r x - r 



+ V xc {r), (18) 



while the nonlocal term generalizes the usual Kleinman- 
Bylander— form allowing several projectors for a given 
angular momentum component 



When augmentation charges vanish (Q nm — 0) the above 
formulas reduce to the standard norm-conserving formu- 
lation. 

In order to generalize Eq. © to the case of US PPs, 
one needs to compute first order perturbation theory in 
presence of overlap matrix, S, as below 
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and AV^ is given by Eq. ©■ 

The derivative of the Kohn-Sham potential is given 
in Ref. 0] and for gradient-corrected functionals in 
Ref. [13. 



III. RESULTS FOR NIOBIUM UNDER 
PRESSURE 

A. Technical details 

The calculations of the ground-state electronic and 
vibrational properties of Nb were performed using the 
Local-Density Approximation (LDA) and an ultrasoft 
pseudopotential. A kinetic energy cut-off of 45 Ry 
(270 Ry) was chosen for the expansion into plane waves 
of the wavefunctions (density). Such high cut-offs were 
necessary to obtain accurate values for some "strategic" 
low-frequency phonons, located mostly near the T-point. 
In fact, even small errors in this region of the spectrum 
lead to large relative errors in the estimate of the a 2 F 
function and of A. 

The integration over the Brillouin zone (BZ) requires 
special techniques to account for the Fermi surface. We 
used the broadening technique proposed in Ref. with 
a smearing parameter of 0.03 Ry (which was tested^ 
to reproduce well the experimental spectra). The grids 
for the electronic BZ integration (k-grid) and for the 
phononic BZ integration (q-grid) have been chosen ac- 
cording to the Monkhorst-Pack schemed 

Details of the numerical quadrature used to evaluate 
the double-delta term appearing in Eq.(2), together with 
convergence tests, are given in the Appendix. 

All calculations were performed using the QUANTUM- 
ESPRESSG^i suite of codes. 
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TABLE I: The lattice constant a (in a.u.), the correspond- 
ing volume ratio V/Vo, the pressure P (in GPa), and bulk 
modulus B (in GPa) calculated for BCC niobium crystal at 
several pressures. The experimental lattice constant— is 6.24 
a.u. and the bulk modulus 23 is 170 GPa. 



B. Structural properties 

Niobium in the body-centered cubic structure was 
studied at eight values of the lattice parameter from 
6.34 a.u. to 5.64 a.u in steps of 0.1 a.u. These lattice pa- 
rameters correspond to pressures ranging from -16.6 GPa 
to 78.4 GPa. The results are reported in TABLE [Q 
The calculated static equilibrium lattice constant (zero- 
point motion and thermal effects not included) is about 
6.14 a.u., slightly underestimating (as it usually happens 
within LDA) the experimental valued of 6.24 a.u. The 
calculated bulk modulus at the theoretical equilibrium 
lattice is 192 GPa, versus a room-temperature exper- 
imental valued of 170 GPa and a calculated value of 
162 GPa at the experimental lattice parameter of 6.24 
a.u. The calculated bulk modulus is very sensitive to 
the volume: it varies by more than a factor three in the 
considered range of pressures. 



C. Band structure and Fermi surface 

The evolution of the band structure as a function of 
pressure is presented in FIG. ^ showing no qualitative 
change in the electronic states at the Fermi surface up to 
about 38.8 GPa. From 56.7 GPa to 78.4 GPa, we observe 
some changes along the T-H and T-N lines. 

A 3D picture of the Fermi surface at ambient pres- 
sure and 56.7 GPa is drawn in FIG. [21 The lower en- 
ergy band forms the octahedron centered at the T-point. 
With increasing pressure, this octahedron shrinks and 
it becomes surrounded by six little ellipsoids when the 
previously described band-structure changes on the Y- 
H line appear. Around N-point, the Fermi surface forms 
ellipsoids that are disconnected at lower pressure and be- 
comes connected by necks to the four neighboring ellip- 
soids above 56.7 GPa. In addition, a complicated open 
sheet structure, often referred as "jungle gym" , extends 
from r to the H points in the BZ. 
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FIG. 1: The band structures of niobium at several pressures 
(in GPa). 



FIG. 2: Fermi surfaces from two bands (the lower energy 
band in the left panels and the higher energy band in the 
right panels) at ambient pressure (top panels) and at 56.7 
GPa (bottom panels) . Pictures obtained with the XCrySDen 
package!— 



Our results are in good agreement with the detailed 
studies of Anderson et alX for the band structure and of 
Ref. 5] for the Fermi surface. 



D. Phonon frequencies and linewidths 

The phonon spectra are presented in FIG. 01 We ob- 
serve an overall increase in phonon frequencies with pres- 
sure, especially at H, P and N high-symmetry points. 
Close to the r-point along the T-H symmetry line, very 
low frequency phonon modes with an anomalous pressure 
dependence can be observed. At both the experimental 
and the calculated equilibrium lattice constants, the two 
transverse modes (Tl and T2 in the following) display 
a very flat dispersion close to r-point. At variance with 
all other modes in the BZ, the Tl and T2 modes along 
T-H line soften with pressure (between «10 to «55 GPa 
in our calculations). This anomalous, non-monotonic, 
dispersion relation eventually disappears and become a 
smooth curve at the highest pressure we have considered. 

The phonon linewidths 7 q „ have a weak dependence 
upon pressure (see FIG.0J, with two important excep- 
tions: i) at low pressure, large variations in phonon 
linewidth occur for the T2 and L modes near the N high- 
symmetry point, and ii) at high pressure (above «55 
GPa), a significant reduction in linewidth is observed in 
many parts of the BZ, especially along the T-H direction. 
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TABLE II: The parameters of niobium under pressure, P (in 
GPa): the electronic density of states at the Fermi surface 
N(ef) (states per spin and per Ry), the electron- phonon cou- 
pling constant A, and the experimental critical temperature 
T c Mp from Ref. Q. 



E. Eliashberg function and electron-phonon 
coupling constant 

By integrating the calculated phonon linewidths and 
frequencies we can obtain Eliashberg a 2 F(ui) function, 
Eq. jnj), that we present in FIG. El for the three phonon 
branches separately. As a general feature, all peaks in 
a 2 F(uj) move to higher frequency with pressure, as ex- 
pected from the global positive frequency shift with pres- 
sure visible in FIG. For the Tl and T2 modes the 
height of the main peak decreases with pressure, while 
for the L mode the maximal height of the peak appears 
at the equilibrium lattice constant. 

Our calculated A and the density of states at the stud- 
ied pressures are presented in TABLED The pressure 
behavior of the electron-phonon coupling constant is sim- 
ilar to what could be expected from the experimental 
pressure dependence for T c : A shows a positive jump at 
low pressure and decreases significantly at high pressure. 

A good review of theoretical works on the electron- 
phonon coupling in Nb at ambient pressure is given by 
Solanki et al. in Ref. |2^, where the reported values of 
the electron-phonon coupling constant range from 0.59, 
obtained 4 from augmented plane- wave method (APW), 
to 1.52, obtained^ later from the same method. More 
recently, SavrasoM^i obtained a value of 1.26, Bauer et 
aim* obtained a value of 1.33, while our results for A at 
ambient pressure is 1.41. 

We notice that the main reason for the spread in the 
reported values of A is the difference in calculated values 
for N(ef), entering the definition of A in the denomina- 
tor (Eqs. ©-0). In fact a large value of 14.1 (states 
per spin and per Ry) for the DOS and a small value of 
0.59 for A have been reported in Ref. [4(, while a value 
of 8.89 for the DOS and A=1.52 have been reported in 
Ref. [2(| . Consistent with this trend, our calculated DOS 
of 10.12 at ambient pressure is somewhat lower than the 
DOS of 10.21 obtained by Savrasovj2£, while other au- 
thors obtained ~ 11 (Ref. ;2]), 11.77 (Ref. |U), and 9.84 
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FIG. 4: Phonon linewidths (in THz) in niobium at several pressures (in GPa). The three modes of 7 q „ correspond to the three 
phonon modes shown in FIG. [3] the scale on the vertical axis is the same for all three modes. 
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FIG. 5: Eliashberg function, a 2 F , in niobium at several pressures (in GPa); the scale on the horizontal axis is the same, within 
a given mode, for all pressures. 
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(Ref. H3). 

The experimental values of A obtained from the elec- 
tronic tunneling spectroscopy^ 3 ^ are 1.04 and 1.22, 
while Haas- van Alphen data^S yield A = 1.33. 



IV. DISCUSSION 

After examination of the phonons and the electron- 
phonon coupling, we notice that the decrease of A at 
high pressure can be easily related to the decrease of the 
a 2 F(u>) peak and to its shift toward higher frequencies 
for all modes. 

The origin of the increase of A at low pressure between 
~0 GPa and «10 GPa is instead more difficult to trace. 
We found out that it is mainly determined by the anoma- 
lous dispersion of the Tl and T2 modes close to T-point 
(see FIG. |3J in the frequency region below 1 THz, that de- 
termines a low-frequency peak in the Eliashberg function 
of the Tl mode above 10 GPa. This peak is instead ab- 
sent at ambient pressure and reappears for an expanded 
lattice only at about -16 GPa. 

It has been shown2i for the Eliashberg model that the 
contribution to T c from acoustic modes close to the T- 
point vanishes. For the low-frequency modes which are 
associated to Kohn anomalies, however, we can expect 
important contribution to T c because phonon softening 
may occur with a finite phonon linewidth. Therefore, 
the region near T can give a very large contribution to 
the electron-phonon coupling in niobium for all studied 
pressures (see TABLE [nj. 

Let us now consider the band structure. In order to 
explain low-pressure anomalies in T c , Struzhkin et alX 
proposed the existence of necks between the ellipsoids 
around N and the "jungle gym" open sheet extending 
from r to H along the T-S line at a pressure below 
5 GPa, and the disappearance of these features at the 
higher pressures. Our calculations do not support this 
suggestion. The detailed analysis of the Fermi surface 
reported by Ostanin et aim- gives results close to ours. 
Previous theoretical 5 i 25 i 26 i 34 i 35 i 36 and experimental 5 ! 37 ! 38 
investigations of the Fermi surface also did not detect any 
changes at low pressure. 

One can, therefore, connect the high-pressure decrease 
in electron-phonon coupling to changes in the band struc- 
ture, while the origin of the low-pressure anomaly re- 
mains unclear. Therefore, we need a different tool in 
order to detect tiny features of the Fermi surface. 

We propose to look closer at Fermi surface nesting by 
plotting the dispersion of the nesting factor: 

X q = ^5(e k -£ F )5(£ k+q -£ F ). (22) 

k 

Large nesting factors correspond to large regions of the 
Fermi surface being connected by the nesting vector, q, 
and are expected to correspond to large electron-phonon 
couplings. 
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FIG. 6: The nesting factor, X q , of the Fermi surface in nio- 
bium at eight studied pressures (in GPa), for selected high- 
symmetry lines (top panel) and high-symmetry points (bot- 
tom panel). 



FIG. 7: The isosurface of the nesting factor in the whole 
BZ, Xq = 1.0 (arbitrary units), for niobium at pressures of 
-0.6 GPa (left panel) and 10.0 GPa (right panel). Pictures 
obtained with the XCrySDen packagei— 



In FIG. El the factor AT q is reported as a function of 
pressure. This quantity has been computed numerically 
using the smearing technique with the broadening of 0.03 
Ry. The maximal nesting takes place at the T-point with 
much smaller maxima at the high-symmetry points H, 
P and N, and in the middle of the lines: r-H, H-P, P- 
r and T-N. The nesting factor decreases monotonically 
with increasing pressure in the whole BZ except around 
ambient pressure (0-10 GPa). At the equilibrium lattice 
constant, a large damping of the nesting factor moves 
the whole nesting curve below its value at 38.8 GPa and 
makes it even similar to the curve drawn for the pressure 
of 56.7 GPa. 

The damping of the nesting factor close to ambient 
pressure occurs in the whole BZ, as we can see in FIG.[7| 
in comparison to nesting for a pressure of 10.0 GPa. This 
damping explains the jump in the total electron-phonon 
coupling constant which happens below and above the 
ambient pressure (see TABLE ITU. 
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V. SUMMARY 



We investigated the origin of the two discontinuities 
of the superconducting critical temperature, observed in 
Niobum at low pressure, about 5 GPa, and at high pres- 
sure, about 60 GPai For this purpose, we developed 
computational tools for the accurate calculation of the 
electron-phonon coupling. 

We find that the anomalous behavior of T c in Nb un- 
der pressure originates in Kohn anomalies close to the 
T-point in the BZ, so the measured discontinuities are 
caused by low- frequency phonons. In agreement with 
previous authors; 2 ' 3 ! 4 ! 5 we find that the high-pressure dis- 
continuity of T c is associated to a visible change in the 
band structure. As for the low-pressure discontinuity, 
we notice that such anomaly shows up as a general de- 
crease of the nesting factor without any visible change in 
the shape of the Fermi surface. Such a decrease is uni- 
form in the whole BZ, explaining why previous calcula- 
tions, Refs. @Q, did not detect any anomaly in the elec- 
tronic structure of Nb near ambient pressure. The total 
electron-phonon coupling constant varies with pressure 
as expected from the measured critical temperatures. 

In conclusion, both discontinuities of T c in niobium, 
at low and high pressures, can be reproduced when 
the electron-phonon spectral function is calculated accu- 
rately, and the anomalies can be explained by the closer 
look into the details of the Fermi-surface nesting and the 
band structure. 
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APPENDIX A: IMPLEMENTATION DETAILS 
AND TEST CALCULATIONS 

The electron-phonon coupling constant A and the 
spectral function a 2 F(ui) are defined by a double-delta 
integration on the Fermi surface (Eqs. J2J) and 
The accurate calculation of these integrands requires a 
very dense sampling in both the electronic (k) and the 
phononic (q) grids. One can use either the broaden- 
ing technique^ or the tetrahedron method4fi We choose 
the former to perform the quadrature on the Fermi sur- 
face, and the latter to evaluate the electron-phonon and 
phonon densities of states as functions of the vibrational 
frequencies. In the broadening scheme, a finite energy 
width is attributed to each state. For any function / 
which has to be integrated with the double-delta, one 




broadening (o") broadening (a) 

FIG. 8: The density of states and the double-delta integrand 
on the Fermi surface of Nb. 

can use the formula 
I = J dk J dq /(k, q) <5(e k - e F )£(e k+q - e F ) 

^M- 1 "^)- (A1) 

where a is the broadening, and N q the number of k- 
and q-points, Qbz the volume of the BZ. For infinitely 
dense grids of k- and q-points, convergence of the in- 
tegrand is achieved when a approaches zero. For finite 
grids, however, one has to find a range of a values yield- 
ing close results for the integrands at different grids. In 
metals like Nb the presence of Kohn anomalies in the 
phonon spectra sets additional requirements for the ac- 
curacy. Thus, it is expected that the aforementioned in- 
tegrands have to be calculated at very dense k-point and 
q-point grids. 

Since the electron-phonon coupling matrix elements 
are smooth functions of k and q, we resort to an interpo- 
lation procedure. For a chosen k- and q-vector grid, we 
calculate matrix elements ffk+q'k™' defined as in Eq. (jSJ 
but with respect to the displacement of a single atom s 
along cartesian component a. For each q-vector we make 
a linear interpolation in k-space of the matrix elements 
to a denser k-vector grid. We then perform the integra- 
tion, using the denser grid and gaussian broadening as in 
Eq. HAljl. of auxiliary phonon linewidths 7^(q), defined 
as: 

— aft { \ \ ^ \ V qsa,mn\* cis'B.mn 

7 a »'(q) = 2^2j5k+q,k ) 5k+q,k 
inn k 

xS(e k+ ^ m - e F )5{e^ n - e F ). (A2) 
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FIG. 9: The phonon linewidth 7 q „ for two selected q-vectors. 
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FIG. 10: The Eliashberg function of Nb for fixed broadening 
(7=0.03 Ry and different q-grid (top panel) and the Eliashberg 
function for fixed q-grid and different broadening a (bottom 
panel) . 
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FIG. 11: The total electron-phonon coupling A of Nb for dif- 
ferent q-grids as a function of broadening a. 
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FIG. 12: The Eliashberg function for Nb at calculated equi- 
librium lattice constant (6.14 a.u.). 



These are related to the 7 q „ of Eq. J^J through the rela- 
tion 

7qv = 27TUV ^^«T^(q)«^. (A3) 

sa s'/3 

Symmetry is exploited to reduce the number of k- and 
q-points used in the calculation. A convenient way to 
achieve such a goal is to perform symmetrization. Let 
us denote with T a " the symmetry operators of the small 
group of q (i.e. the subgroup of crystal symmetry that 
leaves q unchanged) . We restrict summation on k-points 
to the irreducible BZ calculated with respect to the small 
group of q. Then symmetrization is performed, sepa- 
rately for each q-vector, as follows: 



7^(q) = 



J2 T-'T^'rZns^) 

a'j3< 

xexp (iq • (r s - f T(s) - fy + r T( y))) (A4) 



In the above formula, atom s with atomic position r s 
transform into atom T(s) with atomic position tW s ) af- 
ter application of operation T. The symmetrized matrix 
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at each q is subsequently rotated, using the remaining 
crystal symmetries that are not in the small group of q, 
and the symmetrized matrices at all q-vectors in the star 
of q are thus obtained with minimal computational ef- 
fort. The same procedure can be applied to dynamical 
matrices, Eq. (fTTl . 

Once the electron-phonon coupling matrix of Eq. (|A2|) 
are calculated on a q- vector grid, it is possible to perform 
Fourier interpolation and to interpolate to a finer grid. In 
this way, the integration in the q space needed to calcu- 
late a 2 F(uj), Eq. (J2J, and A, Eq. ©, can be accurately 
performed with a reasonable computational effort. 

Let us turn now to some numerical experiments. 
FIG. |H1 shows the density of states for Nb (left panel) 
and the double-delta integrand of a constant function 
(right panel) at the Fermi surface for a selected q-vector, 
as a function of the broadening a and of the Monkhorst- 
Pack 20 k-point grid. For large enough a the results for 
different grids - except the (8,8,8) grid which is too coarse 
- converge to the same values, which however depend on 
a. For small enough a different grids yield different re- 
sults. Since we are interested in the a — > limit, we have 
to choose a grid of an affordable size that yields con- 
verged results for a a as small as possible. A reasonable 
choice is the (64,64,64) grid with cr=0.02 Ry, yielding a 
DOS at the Fermi energy about 10.1 states per spin and 
per Ry. The convergence of the double-delta function is 
a little bit slower than that of a single-delta. 

The phonon linewidths, 7 q „, for two selected phonons 
qv are displayed in FIG. EJ The self-consistent calcula- 



tions were performed i) at the k-grids of (16,16,16) and 
(32,32,32), interpolated to a denser (64,64,64) grid; ii) at 
the (24,24,24) k-grid, interpolated to (72,72,72); iii) at 
the (48,48,48) k-grid, interpolated to (96,96,96). The in- 
tegration weights for the k-space quadrature, i.e. the 
gaussians centered around the single-particle energies, 
were obtained from the accurate self-consistent calcula- 
tion at the corresponding dense grids. As one can see 
in FIG. OH the convergence in k-points is obtained quite 
easily even for the SCF calculation at the grid (16,16,16). 

In order to obtain the spectral function a 2 F(uj) one 
needs to perform the q-space quadrature of the phonon 
linewidths; Eq. J()J). For the phonon and the electron- 
phonon densities of states at given frequency uj, we em- 
ploy the tetrahedron method within the scheme proposed 
by Blochl. 41 FIG. EI| (upper panel) shows the Eliashberg 
function, a 2 F, of Nb for q-grid (8,8,8), without interpo- 
lation and interpolated into denser grids. This has been 
done for a fixed broadening, <j=0.03 Ry. The lower panel 
of the same figure shows the a 2 F function at the q-grid 
of (8,8,8) interpolated to (20,20,20)-point grid for several 
broadenings a. 

In FIG. ^] we report the variation of total electron- 
phonon coupling constant A with the broadening a for 
three q- meshes. 

In FIG. El we present the electron-phonon density of 
states for niobium under ambient pressure. We fit the 
curve of the spectral function with cubic splines^ for a 
finer plot. 
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